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Abstract 

Co-based superalloys in the Co-Al-W system exhibit coherent LI2 Co3(Al,W) 
7' precipitates in an fee Co 7 matrix, analogous to N13AI in Ni-based sys- 
tems. Unlike M3AI however, experimental observations of Co3(Al,W) suggest 
that it is not a stable phase at 1173K. Here, we perform an extensive series of 
density functional theory (DFT) calculations of the 7' Co3(Al,W) phase stabil- 
ity, including point defect energetics and finite-temperature contributions. We 
first confirm and extend previous DFT calculations of the metastability of LI2 
Co3(Alo. 5 Wo. 5 ) 7' at OK with respect to HCP Co, B2 CoAl, and DO19 Co 3 W 
using the special quasi-random structure (SQS) approach to describe the Al/W 
solid solution, employing several exchange/correlation functionals, structures 
with varying degrees of disorder, and newly developed larger SQS. We expand 
the validity of this conclusion by considering the formation of antisite and va- 
cancy point defects, predicting defect formation energies similar in magnitude 
to N13AI. However, in contrast to the N13AI system, we find that substituting 
Co on Al sites is thermodynamically favorable at OK, consistent with exper- 
imental observation of Co excess and Al deficiency in 7' with respect to the 
Co3(Alo.5Wo.5) composition. Lastly, wc consider vibrational, electronic, and 
magnetic contributions to the free energy, finding that they promote the sta- 
bility of 7', making the phase thermodynamically competitive with the convex 
hull at elevated temperature. Surprisingly, this is due to the relatively small 
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finite-temperature contributions of one of the 7' decomposition products, B2 
CoAl, effectively destabilizing the Co, CoAl, and C03W three phase mixture, 
thus stabilizing the 7' phase. 

Keywords: cobalt-base superalloys, point defects, phase stability, density 
functional theory, finite-temperature free energy 
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1. Introduction 



Ni-based superalloys have greatly advanced gas power and jet engine turbine 
efficiencies due to strong high-temperature mechanical performance [1]. This 
strength is commonly attributed to the 7/7' FCC/LI2 microstructure observed 
with Ni and Ni3Al[l]. The discovery of similar 7/7' high-temperature strength- 
ening between FCC Co and LI2 Co3(Al,W)[2] has made Co-based superalloys 
a promising avenue for further improvements in turbine blade performance as 
they have the potential to be stronger and melt at higher temperatures than Ni- 
based systcms[3]. Interestingly, C03AI and C03W in the LI2 structure are not 
stable phases, yet the observed LI2 7', consisting of a disordered solid solution 
of Al and W in the B site of the A3B LI2 structure, is readily observed [2]. The 
reported compositions of 7' include Co-10Al-12W (at.%)[2] from field-emission 
electron probe microanalysis, nearly single-phase 7' cast from Co-10Al-12W and 
Co-10Al-13W[4] and from Co-12Al-llW and Co-10Al-llW[5], and Co-lOAl- 
12. 5W from atom probe tomography [6]. As such, the "ideal" 7' composition is 
often assumed to be equal amounts of Al and W mixing in the B sublattice, 
Co3(Al .5Wo.5), although observed compositions tend to be slightly Co-rich and 
Al-poor relative to this ideal. 

Understanding the thermodynamic stability of the 7/7' two-phase system 
is vitally important for the long-term performance of a Co-based supcralloy. 
The compositional width of the two-phase equilibrium in the ternary diagram 
(i.e. the degree to which the system can tolerate deviations in composition in 
either 7 or 7' before additional phases appear) is fairly narrow[2]. As such, the 
optimization of high-temperature stability and mechanical properties with addi- 
tional alloying elements has been aggressively investigated [3]. However, there is 
evidence that 7' in the Co-Al-W system itself may not be a thermodynamically 
stable phase. Diffusion couples between Co-27Al and Co-15W (at%) initially 
form Co 3 (Al,W) 7' along with FCC Co 7, B2 CoAl, and DO19 Co 3 W) after 
annealing at 1173K for 500 hours[7]. Annealing for 2000 hours, however, in- 
creases the amount of 7, CoAl, and C03W at the cost of 7', suggesting that 7' is 
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metastable with respect to 7, CoAl, and C03W at 1173K. Similar results were 
observed in bulk Co-10A1-7.5W (at%) samples, annealing at 1173K for 2000 
hours [8]. 

Previous studies of Co3(Al,W) 7' with first-principles density functional the- 
ory (DFT) calculations have explored the phase stability, predicting a metastable 
7', with the disordered (^(Alo.sWo.s) LI2 solid solution approximated by a 32 
atom special quasi-random structure (SQS)[9, 10]. LI2 (^(Alo.sWo.s) was 
found to be metastable by 70 meV/atom with respect to the three phase mix- 
ture of FCC Co, B2 CoAl, and DO19 C03WP, 10]. A negative formation energy 
of an ordered LI2 supercell with respect to the elements was also predicted using 
WIEN2K[11]. The effect of additional alloying elements has been investigated 
with DFT as well[10, 12-14]. The DFT predicted single-crystal elastic constants 
of Co3(Al,W) 7' are generally in good agreement with experiments [9, 11, 15]. 

Although the random solid solution of LI2 Co3(Alo.sWo.5) has been pre- 
dicted to be metastable at OK, there are other phenomena which can affect the 
stability of 7'. For instance, ordered structures will be more stable at OK than 
the corresponding random solid solution at the same composition. However it 
is currently unknown how large this ordering energy is for LI2 Co3(Alo.sWo.5). 
Also, the composition of 7' is Co-rich and Al deficient with respect to the "ideal" 
Co3(Al .5Wo.5)[2, 4-6]. Defects altering the composition must be present, such 
as Co vacancies or Al antisite defects, the energetics of which have not been cal- 
culated. Furthermore, finite-temperature phenomena such as vibrational effects 
must be considered to predict whether the OK metastability extends to higher 
temperatures as 7' is observed in samples annealed at high temperature, on the 
order of 1200K. 

As the stability of 7' is such a critical property, we revisit Co3(Alo.sWo.5) 7' 
with DFT, exploring the above-mentioned phenomena which may affect the pre- 
dicted metastability of the phase. We first confirm previous calculations of the 
metastability of Ll 2 Co 3 (Alo. 5 Wo.5) 7' at OK with respect to HCP Co, B2 CoAl, 
and DO19 C03W with the SQS approach, employing several exchange/correlation 
functionals and structures with varying degrees of disorder, including a newly 
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developed larger SQS. All calculations confirm OK LI2 Co3(Alo.sWo.5 metasta- 
bility. The effect of point defects on the stability is then studied, with the 
formation energies of antisitc and vacancy-type defects calculated. We find sim- 
ilar defect formation energies as those for N13AI, except that Co substitutions 
on Al sites are thermodynamically favored. Lastly, finite-tempcraturc contri- 
butions to the free energy of 7' are investigated, including vibrational, thermal 
electronic, and magnetic effects, all of which appear to promote 7' stability, 
primarily due to unique finite-temperature properties of B2 CoAl. 

2. Methodology 

Density functional theory calculations were performed with the Vienna Ab- 
initio Simulation Package (VASP)[16, 17] using the included projector aug- 
mented wave potentials [18]. The exchange and correlation functional was de- 
scribed with the generalized gradient approximation of Perdew, Burke, and 
Ernzerhof[19] (GGA-PBE). The 4s and 3d electrons of Co, the 3s and 3p elec- 
trons of Al, and the 5p, 6s, and 5d electrons of W were treated as valence, unless 
otherwise stated. Crystal structures were relaxed with respect to all degrees of 
freedom with an energy cutoff of 350 eV and gamma-centered k-point meshes 
of approximately 6000 k-points per reciprocal atom. All calculations are spin 
polarized. 

The 7' LI2 Co3(Al,W) disordered solid solution is approximated by special 
quasi-random structures [20], ordered supercclls that reproduce the correlation 
functions of the disordered solution on average over all sites. SQSs have been 
constructed for a variety of common lattice types, including FCC[21], BCC[22], 
HCP[23], and even the situation in the current paper, the B-sublattice of the 
A3B Ll2[9]. Unless otherwise stated, we employ the previously published [24] 
LI2 A3(B,C) mixing SQSs, consisting of 8 mixing atoms (32 atoms total) for 
A3(Bo.5Co.5) and 16 mixing atoms (64 atoms total) for A3(Bo. 25^.75). In the 
current work, we identify SQSs by the number of mixing atoms, for instance 
SQS-8 contains 8 mixing atoms. We also generated an SQS-16 for A 3 (B .5C .5) 
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to better approximate the disordered solution and test the energetic convergence 
of the SQS with respect to supercell size, discussed in detail in the Appendix. 
The OK formation energy of a Co^Al^Wz compound, AE, is calculated by: 

AE(Co x Al y W z ) = E(Co x Al y W z ) - ^ fi m (1) 

m 

where E(Co x Al y W z ) is the DFT total ground state energy of Coa;Al.yW z and 
/i m are the chemical potentials for Co, Al, and W, dictated by the choice of 
reference states. For the formation energy with respect to the elements, AEf, 
the reference states are taken to be HCP Co, FCC Al, and BCC W. For the 
LI2 mixing energy, AE m j x , the reference states are taken to be C03AI and 
C03W in the Ll 2 structure. Of particular interest will be the stability of 7' 
LI2 Co3(Alo.5Wo.5), AEstab, which we define by Equation 1 with respect to the 
set of potentially stable binary and unary structures at this composition: HCP 
Co, DO19 C03W, and B2 CoAl. If AEstab is positive, then 7' is not stable. An 
equivalent relation can be constructed for finite-temperature contributions to 
the free energy of stability of 7' Ll 2 Co3(Al .5Wo.5), AF J gtab : 

AF J stab [T] = F j (Co 3 (Al .6Wo.6))IT] - £ /i™ [T] (2) 

m 

where Fj(Co3(Alo.5Wo.s)) is the contribution from physical effect j to the free 
energy of LI2 Co3(Alo.5Wo.5) and /j,™ is the contribution to the chemical poten- 
tial due to physical effect j . If AF J gtab is negative then j has a stabilizing effect 
on 7'. In this work, the vibrational and thermal electronic contributions to the 
free energy are calculated. The vibrational contribution to the free energy is 
predicted by frozen phonon calculations with the supercell approach in the har- 
monic approximation[25, 26]. Supercell sizes are given in Table 4. The thermal 
electronic free energy is predicted with the Alloy Theoretic Automated Toolkit 
(ATAT)[27]. 

Antisite and vacancy point defects are considered in this work, designated 
by Xy, where X is the species on the Y site of the pristine LI2. □ indicates a 
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vacancy. The defect formation energies, AEDefect(XY), of antisite defects and 
vacancies in the 32 atom SQS-8 structure, C024AI4W4, are calculated by: 

AE Dcfcct = E(Dcfcct) - E(Pristine) - V fi m (3) 

m 

where E(Defect) is the energy of the supercell containing the point defect, 
E (Pristine) is the energy of the perfect C024AI4W4 SQS-8 structure. The chem- 
ical potentials are taken to be HCP Co, B2 CoAl, and DO19 C03W. 

3. Results and Discussion 

3.1. SQS OK Stability 

The formation energy with respect to the elements, AEf, and the Ll 2 mix- 
ing energy, AE m ; x , of the LI2 Co3(Ali_ x W x ) disordered solid solution from 
the SQS approach are shown in Figure 1. Also shown are the energies of a 
series of Al/W ordered LI2 supercells (discussed later) and the predicted ther- 
modynamically stable set of phases at these compositions (the so-called con- 
vex hull), consisting of the HCP Co+B2 C0AI+DO19 C03W three-phase mix- 
ture. Although the Co3(Ali_ x W x ) SQSs are thermodynamically stable with 
respect to the elements, the solid solution is not stable with respect to the other 
phases in the Co-Al-W system at OK. The formation energy of Co3(Alo.sWo.5) 
SQS-8, -125 meV/atom, lies 66 meV/atom (AEstab) above the convex hull. 
In other words, Co3(Alo.sWo.5) LI2 is predicted to transform to HCP Co+B2 
C0AI+DO19 C03W at OK. The SQS results are nearly identical to those of previ- 
ous work [10, 24]. To put the magnitude of this AEstab in context, the configu- 
rational free energy at 1200K due to ideal Al/W mixing is only -18 meV/atom, 
given by — HbT-t (2 In o + \ m \ )- Configurational entropy alone will not make 
Co 3 (Al .sWo.s) more stable than HCP Co+B2 C0AI+DO19 Co 3 W at 1200K. 
1200K will be used as the given temperature for all free energy contributions 
as this is approximately the temperature where 7' was first observed [2] and 
the annealing temperature for the experiments where 7' was observed to be 
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metastable[7, 8]. 

[Figure 1 about here.] 

DFT's ability to predict the stability of 7' can be estimated by comparing to 
experimental formation energies of binary compounds. However, the only exper- 
imentally measured AEf in the Co-Al and Co-W binary systems is for B2 CoAl, 
determined calorimetrically to be -561±30[28, 29] and -694[30] meV/atom. Since 
the current DFT predicted value (-606 mev/atom) lies between the experiments 
and the range between the experiments is so large (130 meV/atom or 20% of 
the DFT value), an assessment of DFT's accuracy in predicting stability in the 
Co-Al-W system is difficult from these experimental values. 

To explore the possibility that the metastability of the LI2 Co3(Alo.5Wo.5) is 
an artifact of some error in the choice of GGA-PBE for the exchange and corre- 
lation functional, AEstab for SQS-8 was calculated with different exchange and 
correlation functionals: GGA-PW91[31] and the local density approximation 
(LDA) of Cepcrlcy and Alder [32] as parameterized by Perdew and Zungcr[33]. 
We also tested a W potential with fewer valence electrons (the 5p electrons are 
placed in the core) with GGA-PBE. The resulting stabilities are given in Table 1. 
The GGA calculations give similar values, varying by only 3 meV/atom. Switch- 
ing from GGA to LDA, on the other hand, has a much larger effect, reducing 
AEstab by 19 meV/atom. However, regardless of DFT settings, Co 3 (Al .5Wo.5) 
SQS-8 is still predicted to be mctastable, suggesting that the prediction of 
metastability is independent of the exchange and correlation functional. 

[Table 1 about here.] 

A key assumption in the current calculations is that the disordered LI2 
Co3(Alo.5Wo.5) solid solution is adequately represented by the SQS-8 structure. 
To test this assumption, AEg ta b for structures of varying degrees of disorder 
were predicted. We have included calculations of the smallest ordered structure 
at the Co3(Alo.5Wo.5) composition, a 2x1x1 LI2 supercell containing alternat- 
ing planes of Al and W (henceforth referred to as the 211 structure), and twenty 



8 



other ordered supercells of varying sizes and Al/W arrangements. We also con- 
structed an SQS twice as large as SQS-8, SQS-16, with correlation functions 
much closer to those of a random solid solution (see Appendix). The result- 
ing AEgtab values, shown in Table 2, range by only 18 meV/atom, suggesting 
that the effect on the total energy due the chemical arrangement of Al and W 
on the B-sublattice of the A3B LI2 is somewhat screened by the Co on the 
A-sublattice. None of the structures are predicted to be stable relative to the 
convex hull, confirming the prediction of 7' metastability at OK is not an ar- 
tifact of the degree of order in the system. The smallest structure is also the 
most energetically favorable, 5 meV/atom more stable than SQS-8. This result 
is consistent with the expectation that ordered compounds are energetically fa- 
vored over disordered solid solutions at OK. The temperature above which the 
Al/W sublatticc will energetically prefer the disordered state over the ordered 
state, from thermal energy alone (fcsT), is 58K (i.e. above this temperature 
the 5 meV/atom difference between the ordered and disordered Co3(Alo.5Wo.s) 
will be overcome by the thermal energy), suggesting that the Al/W sublatticc 
is indeed disordered at the high temperatures 7' is observed. 

[Table 2 about here.] 

3.2. Defects 

As previously stated, relative to the current SQS calculation at Co- 12.5 Al- 
12. 5W, the observed compositions of 7' are typically Co-rich and Al-poor[2, 4-6], 
suggesting defects are present to deviate from the "ideal" Co3(Alo.sWo.5 com- 
position discussed earlier. We investigate deviation from this ideal composition 
by way of antisite defects and vacancies in the SQS-8. These defects have the 
potential to further energetically stabilize the Co3(Alo.5Wo.5) SQS while also 
introducing additional configurational entropy. 

Table 3 summarizes the defect formation energies, energetic stabilities, and 
the configurational free energy (including both defect and Al/W mixing) at 
1200K. A range of energies are given in Table 3 since the SQS-8 contains Al, 
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W, and Co sites of varying local chemistry. All possible defect sites were calcu- 
lated in the current work (e.g. all 24 Co atoms in the SQS-8 were individually 
replaced with Al for the Coai defect). Antisitc defects include replacing Al and 
W with Co, Coai and Cow, respectively, and replacing Co with Al and W, Alc 
and Wco, respectively Notably, Coai has a negative defect formation energy, 
whereas all the other defects are positive. Vacancies of Al, Co, and W (Dai, 
□co, and Dw; respectively) were also calculated. All three vacancy types are 
predicted to be thermodynamically unfavorable, with defect formation energies 
larger than 1000 meV/atom. The formation energies of both the antisite and 
vacancy point defects (except for Coai) are of similar magnitude to their equiv- 
alent in Ni3Al[34]. Also note that the configurational free energy from these 
defects contributes to the free energy. As an example, for both Al/W mixing 
and a single Deo defect in SQS-8 the configurational free energy is calculated 
by -fc B T§(^ln^ + || In ||) - fc B T±(±ln± + ±ln±). With point defects, 
the configurational contribution to the free energy at 1200K increases from - 
18 meV/atom for only ideal Al/W mixing to, at most, -31 meV/atom. Single 
site defect calculations were performed on the SQS-16 structure as well to test 
whether the SQS-8 supcrcell size is sufficiently large to isolate the point defects 
and minimize artificial defect-defect interactions due to the periodic boundary 
conditions employed in DFT. The SQS-16 defect formation energies (shown in 
parentheses in Table 3) are similar to those of SQS-8, indicating the current 
SQS-8 defect formation energies are converged with respect to supercell size. 

[Table 3 about here.] 

The negative defect formation energy for Coai indicates that replacing Al 
with Co stabilizes the SQS-8 energetically with respect to the convex hull at OK, 
perhaps approaching a more stable Co-rich ordered compound. This suggests 
7' may prefer Co-rich compositions from the ideal Co3(Alo.sWo.5), in agree- 
ment with the observed 7' compositions[2, 4-6]. Furthermore, the observed W 
composition is very similar to the ideal 12.5 at.%, consistent with the current 
predictions since Cow is not as competitive as Coai- To explore the limits of 
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Coai favorability, additional Al atoms were substituted with Co (indicated as 
2Coai, 3Coai, and 4Coai for two, three, and four Al replacements, respectively). 
The results are given in Table 3, where the defect formation energies are taken 
with respect to the previous Coai defect. Interestingly, the second Coai replace- 
ment is also favorable, by about the same energy as the first. The third and 
fourth, however, are no longer favorable. Although Coai defects are favorable, 
their effect on the stability of 7' is small, reducing the metastability by less 
than 10 meV/atom for 2Coai- In the A3B LI2 structure, a B-site has six near- 
est neighbor B-sites. Al-sites that correspond to favorable Coai substitutions in 
the SQS-8 happen to be those that are rich in Al in these six nearest neighbor 
B-sitcs, indicating Co-Al LI2 nearest neighbor B-site bonds are favored. 

3.3. Finite- Temperature Effects 

Beyond configurational disorder, other finite-temperature contributions can 
alter the stability of 7' at the elevated temperatures in which it is experimen- 
tally observed, including vibrational, electronic, and magnetic effects. Of these 
phenomena, phonons, the characteristic vibrations of a crystalline lattice, can 
have a large impact on the stability of metal alloys [35-39]. If the vibrational 
properties of 7' are sufficiently different from the competing phases, then there 
is the possibility the vibrational free energy of stability, AFg^, will be negative, 
having a stabilizing effect on 7'. AFg 7 ^ for Co3(Alo.5Wo.s) is predicted from 
DFT by calculating, with the frozen phonon supercell approach, the vibrational 
free energies of HCP CO, B2 CoAl, DO19 Co 3 W, and Ll 2 Co 3 (Alo. 5 Wo.5) for 
the 211 and SQS-8 structures. The vibrational free energies of these structures 
at 1200K are given in Table 4. The 16 atom Co vibrational free energy is less 
than 1% different from the 54 atom value, indicating 54 atoms is sufficiently 
converged for the current work. Interestingly, the difference between the SQS-8 
and 211 ordered vibrational free energies is also very small, suggesting that the 
vibrational properties of the chemical ordering of the B sublattice in A 3 B Ll 2 
are screened by the A sublattice. AFg 7 ^ at 1200K is -27 meV/atom, contribut- 
ing a stabilizing effect to Co3(Alo.5Wo.s). The negative AFg^ is due to the 
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vibrational free energy of CoAl, which is smaller than those of the other struc- 
tures, indicating stiffer bonds than the other structures. A calculation of the 
CoAl vibrational entropy with a larger supcrcell results in a similar value (also 
given in Table 4), suggesting the CoAl vibrational free energy is converged with 
respect to supercell size. 

[Table 4 about here.] 

Another possible contribution to the free energy is the thermal excitation of 
the electrons from their OK ground state below the Fermi level to excited states. 
This contribution must be considered for metals since the electronic density of 
states (DOS) will be non-zero at the Fermi level[40]. The electronic DOS, pre- 
dicted from DFT, for HCP Co, B2 CoAl, D0i 9 Co 3 W, and Ll 2 Co 3 (Alo. 5 Wo. 5 ) 
with the 211, SQS-8, and SQS-16 structures are shown in Figure 2. Electronic 
excitations will affect 7' stability if the 7' DOS near the Fermi level is different 
from those of the competing phases. This difference is calculated by taking the 
weighted difference of the DOSs as in AEstab (effectively a DOS of stability), 
using SQS-8 for Co3(Al .5Wo.5). Also shown in Figure 2, the DOS of stability is 
non-zero near the Fermi level, particularly just below it. This peak in the DOS 
of stability is due to B2 CoAl, which exhibits a deep well in the DOS below the 
Fermi level, in agreement with previous calculations [41]. Calculation of the elec- 
tronic free energies at 1200K, F^ookj are gi ven m Table 5 and confirm the imbal- 
ance between 7' and the other structures. The DOSs for the LI2 Co3(Alo.5Wo.s) 
structures show good agreement with previous calculations[9, 11], and the dif- 
ference in the electronic free energies and DOSs between the three structures 
is minor, suggesting that the arrangement of Al and W atoms has little effect 
on the electronic free energy. The electronic free energy of stability, AFg t ' Qb , 
calculated from the SQS-8 electronic free energy, is -10 meV/atom. In other 
words, the electronic free energy also has a stabilizing effect on Co 3 (Al .5Wo.5). 

[Table 5 about here.] 

[Figure 2 about here.] 
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Accurate prediction of magnetic thermodynamics at elevated temperatures 
with DFT requires sophisticated models [42-48] beyond the scope of the current 
work. However, some general comments can be made from the current OK 
ground state calculations. As with vibrational effects, magnetism can affect 
the stability of 7' if the any of the phases it competes with have sufficiently 
different magnetic properties. Co, C03W, and 7' SQS-8 are predicted to have 
ferromagnetic order on the Co sites, with average magnetic moments of 1.61, 
0.31, 0.74 [1b, respectively. CoAl, however, is diamagnctic in DFT, with no local 
Co moment. This is an experimentally observed phenomenon[41], attributed to 
the Al valence electrons pairing with the Co valence electrons responsible for 
magnetism. Therefore, the magnetic free energy of stability of 7' is most likely 
negative, as CoAl has no magnetic entropy to contribute towards stabilizing 
the HCP Co+B2 C0AI+DO19 C03W three-phase mixture. An upper bound on 
the sensitivity of AEstab due to magnetism can be estimated by calculating the 
stability from DFT total energies without spin moments considered, essentially 
the nonmagnetic 7' OK stability. We predict this value to be 12 meV/atom, 56 
meV/atom lower than the magnetic value, indicating that magnetism has the 
potential to greatly affect to the stability of 7'. 

Taken altogether, the calculated contributions to the free energy of sta- 
bility for 7' Ll 2 Co3(Al .5Wo.5) are predicted to promote the stability of the 
phase. At 1200K, the free energy contribution from the ideal configurational 
entropy of Al/W mixing is -18 meV/atom, -27 meV/atom from vibrations, and 
-10 meV/atom from thermal electronic excitations. These effects reduce the 
metastability of the SQS-8 structure from 66 meV/atom above the convex hull 
to only 11 meV/atom. Additional stabilization from defects and magnetism 
will likely reduce it further, possibly making 7' thermodynamically competitive 
with the other phases in the Co-Al-W system. 
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4. Summary 

The DFT prediction of 7' LI2 (^(Alo.sWo.s) metastability has been ex- 
plored in detail. Various structures with different degrees of disorder and several 
DFT exchange/correlation functionals agree that Ll 2 Co3(Al .5Wo.5) lies above 
the convex hull of HCP Co, B2 CoAl, and DO19 Co 3 W at OK, with the PBE 
SQS-8 prediction as metastablc by 66 meV/atom. Substitutional point defects 
have small positive defect formation energies while vacancies are much more 
unfavorable. However, substituting Al sites with Co is predicted to be thermo- 
dynamically favorable at OK, in agreement with experimental measurements of 
excess Co and Al deficiency in 7' relative to the ideal Co3(Alo.5Wo.s) composi- 
tion. Finite-temperature contributions to the free energy have been predicted as 
well. Lattice vibrations and electronic excitations appear to have a stabilizing 
contribution to the free energy of LI2 Co3(Alo.5Wo.s) due to small free energy 
contributions from B2 CoAl. Magnetism may also promote 7' stability as CoAl 
is nonmagnetic. Summing all the predicted contributions to the free energy 
of 7' reduces the metastability from 66 to 11 meV/atom. 7' may become less 
metastablc, possibly even stable, when considering further finite-temperature 
contributions, such as magnetism. 
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Appendix: Simple Cubic SQS-16 

The disordered solid solution on the B-site in the A3B LI2 of 7' is described 
with special quasi-random structure approach[20]. The B-site of LI2 is equiv- 
alent to a simple cubic lattice. An SQS for this sublattice at A3(Bo.5Co.s) 
composition containing 8 mixing atoms (SQS-8) has previously been developed 
and employed to test the stability of Co3(Al,W) 7' [9]. To assess the ability 
of this SQS-8 to predict the enthalpy of mixing for 7', we generate a simple 
cubic SQS with twice the mixing sites, SQS-16, which will serve as a better 
approximation for the disordered solid solution. 

We begin by enumerating all possible SQS structures of up to 16 mixing 
atoms (360,642 structures) using a lattice enumeration algorithm [49]. The pair, 
triplet, and quadruplet correlation functions, a measure of the degree of chemical 
disorder in a fixed lattice, for all 360,642 structures were calculated with the 
corrdump function of ATAT[27]. Of all possible structures, a single structure, 
detailed in Table 6, has pair correlations identical to that of the random solution 
out to the 14th nearest neighbor shell. The next best structure approximates 
the random solutions only as far as the 6th. In comparison, the published 
SQS-8 deviates from the random solid solution pair correlation function by the 
4th nearest neighbor shell. On this basis, the structure in Table 6 was chosen 
as the SQS-16 used in the DFT calculations in the current work. The pair, 
triplet, and quadruplet correlation functions for the new SQS-16, the published 
SQS-8, and the ideal disordered solid solution are given in Table 7. The new 
SQS-16 performs better than the published SQS-8 at approximating the pair 
and quadruplet correlation functions of the disordered solid solution, whereas 
similar triplet correlation performance is present in the two SQSs. 

[Table 6 about here.] 

[Table 7 about here.] 
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Table 1: Energetic stability of LI2 Cc^Alo.sWo.s) SQS-8 with respect to the 
convex hull, in meV/atom, for various exchange and correlation functional and 
W potentials. 



A Estab 

66 
63 
47 
65 



GGA-PBE 

GGA-PW91 

LDA 

GGA-PBE (W core 5p) 
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Table 2: Energetic stability, in meV/atom, of Ll 2 Cc^Alo.sWo.s) structures 
with varying supercell size and disorder, with the GGA-PBE potential. 



Structure 


# Al/W 


AEstab 


211 


2 


61 


20 Other Ordered 


4-8 


62-79 


SQS-8 


8 


66 


SQS-16 


16 


67 
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Table 3: Energetics of antisite defects and vacancies in Cc^Alo.sWo.s) SQS- 
8 (and SQS-16 in parentheses), including defect formation energy, energetic 
stability with respect to the convex hull, and the configurational free energy 
(including both defect and Al/W mixing) at 1200K. NCoai refers to N additional 
substitutions of Al with Co, and their corresponding AEDcfcct is for the cost of 
placing the additional defect. 





AEDefect 

[meV/defect] 


AEstab 

[meV/atom] 


pCJontig 
r 1200K 

[meV/atom] 


Com 


-93 - 25 (-115) 


64 - 67 


-25 


2Coai 


-91 - 63 


61 - 66 


-27 


3Co A i 


101 - 102 


64 - 64 


-25 


4Co A i 


162 


69 


-18 


Cow 


275 - 423 (285) 


75 - 80 


-25 


Alco 


263 - 825 (409) 


75 - 92 


-31 


Wco 


300 - 1464 (857) 


76 - 112 


-31 


□ai 


1996 - 2001 (2034) 


133 - 133 


-25 


□w 


2255 - 2438 (2312) 


141 - 147 


-25 


□Co 


1476 - 1932 (1637) 


116 - 131 


-31 
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Table 4: Vibrational free energies, in meV/atom, at 1200K predicted by DFT 
with the supercell phonon approach for Co-Al-W structures. 



# Atoms 

HCP Co 16 -451 

54 -449 

B2 CoAl 54 -361 

128 -365 

D0i 9 Co 3 W 64 -443 

Ll a Co 3 (Al . 5 W . 5 ) 

211 54 -451 

SQS-8 32 -453 
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Tabic 5: Electronic free energy at 1200K predicted by DFT, in meV/atom. 



~5E1 

r 1200K 



HCP Co -16 
B2 CoAl -8 
DO19 Co 3 W -12 



Co 3 (Al . 5 Wo. 5 ) 

211 -23 

SQS-8 -22 

SQS-16 -23 
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Table 6: Crystal structure of the simple cubic SQS-16 developed and employed 
in the current work. 











0\ 






Lattice Vectors 


(:!- 


2 1 












-2 3/ 




Element 1 






Element 2 




X 


y 


z 


X 


y 


z 


0.000 


0.000 


0.000 


0.500 


0.625 


0.625 


0.500 


0.125 


0.125 


0.000 


0.750 


0.750 


0.000 


0.250 


0.250 


0.500 


0.875 


0.875 


0.500 


0.375 


0.375 


0.500 


0.500 


0.000 


0.000 


0.500 


0.500 


0.000 


0.625 


0.125 


0.500 


0.750 


0.250 


0.000 


0.875 


0.375 


0.000 


0.125 


0.625 


0.500 


0.000 


0.500 


0.000 


0.375 


0.875 


0.500 


0.250 


0.750 
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Table 7: Simple cubic correlation functions, IIfe iTO [n], of pairs (k = 2), triplets 
(k = 3), and quadruplets (k = 4) at the mth nearest neighbor shell, each with 
n equivalent figures. Correlation functions are provided for the random solid 
solution, the SQS-16, and the previously published SQS-8[9]. 
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1 


n 2 , 4 [i2] 








-1/3 


n 3>3 [3] 











n 4)2 [24] 











n 2 , 5 [i2] 











n 3 , 3 [i2] 











n 4>2 [24] 





-1/12 





n 2)6 [6] 











n 3)3 [i2] 











114,2 [6] 





1/6 


-1/3 


n 2 , 7 [3] 








-1/3 


n 3)4 [24] 











n 4 . 3 [i2] 











n 2 , 7 [i2] 











n 3 , 4 [24] 











n 4 , 3 [i2] 





-1/6 





n 2 , 8 [i2] 











n 3 , 4 [48] 











n 4 , 3 [3] 





1/3 


-1/3 


n 2 , 9 [i2] 











n 3 , 4 [24] 











n 4 , 3 [i2] 





1/6 





n 2)9 [6] 











n 3 . 4 [i2] 











114,3 [6] 








-1/3 


n 2 ,i [i2] 








1 


n 3 . 4 [i2] 











n 4 . 4 [24] 











n 2 ,u[i2] 











n 3 , 5 [48] 











n 4 , 4 [i2] 





1/6 


-1/3 


n 2 ,i 2 [24] 











n 3 , 5 [24] 











n 4 , 4 [24] 





1/6 





n 2 ,i 3 [i2] 





1 


1 


n 3>5 [24] 











n 4>4 [48] 





-1/120 






24 





Figure 1: DFT formation energies with respect to the elements (top) and mixing 
energies with respect to the C03AI and C03W Ll 2 structures (bottom) for the 
Ll 2 Co 3 (Ali_ x W x ) SQSs, the Ll 2 ordered structures given in Table 2, and the 
C03AI-C03W convex hull, consisting of HCP Co+B2 C0AI+DO19 Co 3 W. 



25 



o 

> 
U 



t/3 
O 
Q 



o 
> 

C3 



t/3 

o 

Q 



o 

> 

c/3 



00 

O 
Q 



o 

> 

u 

-f— < 

C3 



C/3 

o 
p 



o 

> 

1) 
3 



on 
O 
Q 




Co 



-10 -5 5 10 

E-E F [eV] 




E-E F [eV] 




5 

E-Ep [eV] 



10 





Stability 







-10 -5 5 10 

E-E F [eV] 

26 

Figure 2: Total electronic densities of states for HCP Co, B2 CoAl, D0i 9 Co 3 W, 
and Ll 2 Co 3 (Alo. 5 Wo.5) with the 211, SQS-8, and SQS-16 structures. The 
difference between the Ll 2 Co3(Alo.5Wo.s) SQS-8 DOS and the other phases is 
shown as a DOS of stability, calculated by taking the weighted difference of the 
DOSs as in AEg ta b- Energies are normalized to the Fermi level (Ep). 



